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Motivated by the unconventional properties and rich phase diagram of Na^CoC^ we consider the 
electronic and magnetic properties of a two-dimensional Hubbard model on an isotropic triangular 
lattice doped with electrons away from half-filling. Dynamical mean-field theory (DMFT) calcula- 
tions predict that for negative inter-site hopping amplitudes (t < 0) and an on-site Coulomb repul- 
sion, U, comparable to the bandwidth, the system displays properties typical of a weakly correlated 
metal. In contrast, for i>0a large enhancement of the effective mass, itinerant ferromagnetism and 
a metallic phase with a Curie- Weiss magnetic susceptibility are found in a broad electron doping 
range. The different behavior encountered is a consequence of the larger non-interacting density of 
states (DOS) at the Fermi level for t > than for t < 0, which effectively enhances the mass and 
the scattering amplitude of the quasiparticles. The shape of the DOS is crucial for the occurrence of 
ferromagnetism as for t > the energy cost of polarizing the system is much smaller than for t < 0. 
Our observation of Nagaoka ferromagnetism is consistent with the A-type antiferromagnetism (i.e. 
ferromagnetic layers stacked antiferromagnetically) observed in neutron scattering experiments on 
Na K Co02. The transport and magnetic properties measured in Na I Co02 are consistent with DMFT 
predictions of a metal close to the Mott insulator and we discuss the role of Na ordering in driv- 
ing the system towards the Mott transition. We propose that "Curie- Weiss metal" phase observed 
in Na :E Co02 is a consequence of the crossover from "bad metal" with incoherent quasiparticles at 
temperatures T > T* and Fermi liquid behavior with enhanced parameters below T* , where T* is 
a low energy coherence scale induced by strong local Coulomb electron correlations. Our analysis 
also shows that the one band Hubbard model on a triangular lattice is not enough to describe the 
unusual properties of Na :E Co02 and is used to identify the simplest relevant model that captures 
the essential physics in Naa;Co02. We propose a model which allows for the Na ordering phenomena 
observed in the system which, we propose, drives the system close to the Mott insulating phase even 
at large dopings. 

PACS numbers: 71.10.Hf, 71.30.+h, 



I. INTRODUCTION 

One of the outstanding problems in quantum many- 
body physics is to understand quasi-two-dimensional sys- 
tems in which both electron-electron interactions and 
frustration effects are strongi New exotic phases are ex- 
pected to result from the interplay of these effects. Re- 
cently a number of materials have been discovered where 
the interplay of geometrical frustration and strong elec- 
tronic correlations lead to the emergence of unconven- 
tional phases. Examples of such materials include in- 
organic materials and organic charge transfer salts such 
as k-(BEDT-TTF) 2 Cu 2 (CN) 3 or /3'-[Pd(dmit) 2 ]X (Refs. 
and explain this nomenclature) . Antiferromagnets 
on triangular lattices may be isotropic as in NiGa 2 S4, n- 
(BEDT-TTF) 2 Cu 2 (CN) 3 , and, the material we will focus 
on, Na x Co0 2 , or anisotropic as in Cs 2 CuCLi, k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]C1. There is also a great deal of inter- 
est in other frustrated geometries such as the Kagome 
and pyrochlore lattices. However, NazCo0 2 has the, 
so far, unique property among strongly correlated tri- 
angular lattice compounds, that it can be doped. In 
k-(ET) 2 Cu 2 (CN) 3 a spin liquid state gives way to su- 
perconductivity when hydrostatic pressure is applied. 4 



The /3'-[Pd(dmit) 2 ]AT series^ of compounds, Cs 2 CuCl4 
(Ref. @) and NiGa 2 S 4 (Ref. 0) may also show spin liq- 
uid behavior. Interestingly, the underlying structures of 
k-(ET) 2 Cu 2 (CN) 3 and the /3'-[Pd(dmit) 2 ]X compounds 
are very similar to that of NiGa 2 S4 and Na x Co0 2 with ei- 
ther the organic molecules, the Ni atoms or the Co atoms 
arranged in a triangular geometry. Understanding the 
physics of these systems in which both frustration and 
correlations are present is clearly a major challenge to 
theory. 

Much attention has been focused on the occurrence of 
superconductivity in Na a; Co0 2 -yH 2 0. However, the non- 
superconducting, non-hydrated counterpart, Na^CoO^ 
displays many interesting magnetic and electronic prop- 
erties. These include a charge ordered insulating state 
in a narrow band of doping around^ x = 0.5, A-type 
antiferromagnetism 9 for x > 0.75 and a metallic phase 
with a Curie- Weiss susceptibility, x(T) ~ 1/(T + 0), 
whose magnitude is much larger than the Pauli suscep- 
tibility expected for a weakly interacting metal^ (the 
Curie- Weiss metal). This phase also has an extremely 
large thermopower (of order fee /e) that is temperature 
and field dependent i& Foo et al. suggested^ that the 
metallic state for x < 0.5 should be regarded as param- 
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agnetic and the state for 0.5 < x < 0.75 be thought of as 
a Curie-Weiss metal. These results should be compared 
with the measurements of Prabhakaran et aim* which 
display a different behavior with the Curie- Weiss behav- 
ior strongest at x = 0.5 (the lowest value of x in their 
sample set) and the susceptibility becoming less tem- 
perature dependent as x is increased. These differences 
may be due to the rather different strength fields used 
in the two different experiments^ a strong dependence 
on the direction of the field or to some extrinsic effect. 
For x > 0.75 A- type antiferromagnetism (in-plane ferro- 
magnetic ordering stacked in alternate directions along 
the c-axis) is observedi^i^iSiii On the triangular lat- 
tice this raises the interesting possibility of in-plane Na- 
gaoka ferromagnetism^ which we will discuss further 
below. There is also direct experimental evidence for 
strong electron-electron correlations which comes from 
unconventional behavior of transport properties. For in- 
stance, a Fermi liquid like resistivity \p(T) = po + AT 2 ] 
is only observed^ for T < 1 K, and the value of the 
Kadowaki- Woods ratio is comparable to ruthenates and 
heavy fermionSfi2i211 The low temperature scale and the 
Kadowaki- Woods ratio are strongly field dependent with 
the system becoming more weakly correlated with in- 
creasing field. 

Interestingly, some of the above properties are rem- 
iniscent of the predictions of the dynamical mean-field 
theory (DMFT) of the Hubbard modeliSi At half filling 
(n = 1) DMFT predicts a Mott metal-insulator transi- 
tion for U ~ W . In the metallic phase close to the Mott 
transition there is a low temperature scale, T*(n), at 
which there is a smooth crossover from a Fermi liquid to 
an incoherent "bad metal" Transport properties 

display a crossover from those expected in a Fermi liquid 
to properties characteristic of incoherent excitations as 
the temperature is increased. Within the "bad metal" 
phase, that is for T > T* , electrons behave as quasi- 
localized moments which leads to a Curie- Weiss behavior 
of the magnetic susceptibilityiSiiSS The above discussion 
suggests that DMFT, which treats the local electronic 
correlations exactly, may capture the relevant physics 
needed to describe the temperature dependence of many 
of the transport and magnetic properties of Na x Co02. 

An important first step in studying the properties of 
Na x Co02 is to consider the two-dimensional triangular 
lattice Hubbard model. This is motivated by the ar- 
rangement of the Co atoms and because correlation ef- 
fects can be singled out and better understood in the 
simplest model. Related models such as the t — J model 
on the triangular lattice have been recently analyzed us- 
ing resonating valence bond (RVB) theory^ The implicit 
assumption in this kind of approach is that many impor- 
tant aspects of more realistic models for Na x Co02 are 
already contained in the simpler triangular lattice model. 
Our present study, however, leads to the important con- 
clusion that a single band Hubbard model on a triangular 
lattice is not enough to describe the magnetic and elec- 
tronic properties of Na2,Co02 and a more realistic model 



is proposed. 

Our main finding, summarized in Fig. ^ is that DMFT 
predicts dramatic changes in the behavior of the Hub- 
bard model on a triangular lattice when the sign of t 
is changed. For t > we find significant effects due to 
strong electronic correlations. We find that in the elec- 
tron doped triangular lattice a Curie- Weiss metal and 
metallic ferromagnetism arise when only local Coulomb 
correlations are taken into account. The Curie- Weiss 
metal occurs for sufficiently large Coulomb repulsion en- 
ergies U >W and t > 0, where W is the bandwidth. In 
this parameter region, the local magnetic susceptibility 
smoothly changes from the Curie form characteristic of 
local moments at high temperature to Fermi liquid be- 
havior at low temperatures, T < T*(n) <C ep, where cf is 
the Fermi energy. In the Curie- Weiss metal the uniform 
susceptibility is strongly enhanced by both band narrow- 
ing due to correlations and the proximity to the ferro- 
magnetic instability. In contrast, the electron doped tri- 
angular Hubbard model with t < and U > W displays 
Pauli paramagnetism weakly enhanced by the Coulomb 
repulsion; a behavior typical of weakly correlated metals. 
Thus our analysis shows how the non-interacting disper- 
sion has a dramatic effect on the electronic and magnetic 
properties of a frustrated Hubbard model. Thus, our 
work is relevant to DMFT applications on real materials 
where similar LDA density of states appear helping to 
understand the physics in the more complicated situa- 
tions. 

The remainder of this paper is organized as follows. In 
Sec. |Ii|we introduce the Hubbard model on a triangu- 
lar model discussed above, and then examine the solu- 
tion of this model on three sites which already contains 
some features of, and gives significant insight too, the 
solution of the triangular lattice in the thermodynamic 
limit. In Sec. EH we introduce the main DMFT equa- 
tions and present arguments which suggest that DMFT, 
which is only exact in infinite dimensions, may capture 
the essential physics of the Hubbard model on a two- 
dimensional triangular lattice due to frustration. In Sec. 
IIVI we analyze the temperature dependence of the local 
and uniform magnetic susceptibilities for both t < and 
t > 0. In Sec. 0we analyze the local correlation effects 
and the electronic properties of the model. We report 
spectral densities and self-energies in both signs of t and 
their relation with the different magnetic susceptibilities 
found in Sec. IIVI In Sec. I VII we show that for t > and 
electron doping we find Nagaoka ferromagnetism. In Sec. 
IVIII we discuss the experimental situation on Na^Co02 
in the light of our results. In Sec. IVIIII we discuss im- 
provements and extensions to both the model and the 
approximation used to study it, we argue that multiple 
band models and particularly the effects of ordering of 
the Na ions are essential to understand the behaviour ob- 
served in Na x Co02. The paper ends in Sec. IIXI in which 
a summary of main results is given and our conclusions 
are drawn. 
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FIG. 1: Phase diagram for the electron doped Hubbard model 
on a triangular lattice obtained from DMFT calculations. 
U/\t\ is the ratio of the on-site Coulomb repulsion to the ab- 
solute value of the hopping matrix element t and n is the av- 
erage number of electrons per site. The sign of t has a leads to 
qualitative changes in the phase diagram. In the t > case, a 
ferromagnetic metal and a metal with Curie- Weiss suscepti- 
bility are found. In contrast, for !<0a paramagnetic metal 
appears throughout the phase diagram. To make a direct 
comparison with Na^CoO^ one should note that x = n — 1. 



II. THE HUBBARD MODEL ON A 
TRIANGULAR LATTICE 

As stated in the introduction, we analyze the magnetic 
and electronic properties of a Hubbard model on a trian- 
gular lattice, 



H = —t 



( C icr C j<? 



B.c. 



(i) 



where fj,, t and U are the chemical potential, the nearest- 
neighbor hopping amplitude and the on-site Coulomb re- 



pulsion energies, respectively, c, 
electron on site i with spin a, 



(t) 



(creates) destroys an 
. ) indicates that the 



sum is over nearest neighbours only and the number op- 
erator is n irT = c\ a c ia . 

For [7 = the above model gives a tight-binding dis- 
persion 



£k = — 2icos(fc x ) — 4icos cos (^fj 



(2) 



with bandwidth W = 9\t\ 



A. Review of previous work 

For the Hubbard model on the square lattice at half 
filling, the ground state is believed to be a Mott insulator 
with Neel order for all values of U/t. The fact that this 
occurs even for arbitrarily small U /t is due to the per- 
fect nesting of the non-interacting Fermi surface at half 
filling. In contrast, the triangular lattice may exhibit 
diverse phases as U/t is varied. Possible ground states 
that have been proposed include the Mott insulator, com- 
mensurate non-collinear antiferromagnetism, incommen- 
surate spin-density wave (both metallic and insulating), 
spin liquid, superconducting, and metallic states. 

In the strong correlation limit, U t, there is a gap in 
the charge excitation spectrum and the ground state is 
a Mott insulator. The spin excitations are described by 
an antiferromagnetic Heisenberg model on the triangular 
lattice with exchange constant J = At 2 /U. Numerical 
calculations from a range of techniques including exact 
diagonalisation, 24 variational quantum Monte CarlOfS&iSSi 
coupled cluster methods, 27 and series expansions' 8 ^ 
suggest that this has a magnetically ordered ground state 
in which neighbouring spins are rotated by 120 degrees 
relative to one another, as in the classical ground state. 
This agrees with the predictions of spin wave theory^ 
However, series expansions calculations suggest that the 
excitation spectrum is qualitatively different from that 
predicted by non-linear spin wave theory»2i 

At half-filling there is no nesting of the non- interacting 
Fermi surface and so for small U/t a metallic state is pos- 
sible. Exact diagonalisation calculations of the Drude 
weight for lattices with 12 sites suggest that there is 
a first-order transition from a metal to a Mott insu- 
lator when U ~ 12.1 i, 3 - mean-field RVB calculations 
in the thermodynamic limit find a remarkably similar 
result with a first order Mott transition occurring at 
U ~ 12.4 tM When U ^ At the effective spin Hamil- 
tonian in Mott insulating phase should include ring ex- 
change terms of order t 3 /U 2 (Ref. 34). On the triangular 
lattice this can lead to a spin liquid ground stated The 
possible realisation of a spin liquid in k-(ET)2Cu2(CN)3 
(Ref. 0) has stimulated the proposal of specific analytical 
theories for the ground state in this regimeiSS^ 7 . In this 
parameter regime a superconducting ground state with 
d + id symmetry, 38 or odd frequency pairing 39 have also 
been proposed [it is therefore interesting to note that odd 
frequency pairing has also been discussed on the basis 
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of specific microscopic calculations for Na^OoO^ • 2/H2O 
(Ref. E3)]. 

Hartree-Fock^i and slave boson mean-field 
calculations— have found that near the metal-insulator 
phase boundary that incommensurate spin-density wave 
states are also stable. Simulated annealing calculations 
of mean-field solutions with large possible unit cells 
found that for U — 2t and U = At that the most stable 
states were metals with charge and spin ordering. 43 

For electron doping away from half-filling and U t, 
doubly occupied states can be projected out and at — J 
model can describe the system. Extensive studies of this 
model have been madei 12 i 44 These studies show that the 
sign of t has a significant effect on the type of ground 
stated RVB calculations 12 ' 44 found that, for t > 0, d+id 
superconductivity was possible for a range of dopings 
and suggested a variety of unusual metallic states. Series 
expansions^ have emphasised the competition between 
singlet and triplet formation on the triangular lattice in 
terms of competition between RVB states and Nagaoka 
ferromagnetism. This work concludes that for hole dop- 
ing the RVB state is favoured by t > and the Nagaoka 
state is favoured by t < 0. Changing the sign oft is equiv- 
alent to changing from electron to hole doping, therefore 
this result is entirely consistent with the results sum- 
marised in Fig. ^ A more detailed comparison between 
our results and those of Koretsune and Ogata^ will be 
given in Sec. IVfl 



B. Relevance to Na x Co02 

For Na x Co02 the LDA bandwidth gives \t\ ~ 
0.1 eVM^ For electron doping the dispersion @ leads 
to a Fermi surface which is hole-like for t < and 
electron-like for t > 0, thus t < is roughly con- 
sistent with the band structure calculations and angle 
resolved photoemission spectroscopy (ARPES) experi- 
ments for Na^CoC^. However, the band structure of 
Na x Co02 is still controversial. ARPES 46 " seems to sug- 
gest that a simple tight binding model with t < on 
the triangular lattice may be sufficient to describe the 
dig bands of Na^OoO^. However, even within a single 
band model hopping matrix elements up to at least 3rd 
nearest neighbors 4 ^ 4 *!*^ need to be included in order to 
correctly reproduce the LDA band structure. Further, 
LDA calculations^^ 8 , predict the existence of six small 
elliptical hole pockets near to the K-points if x < 0.6. ft 
is not clear, at this point, whether or not the elliptical 
pockets actually existi^ 1 *^ In contrast to the a% g bands, 
the LDA dispersion of the e\ g bands is better fitted to a 
t > triangular lattice tight-binding dispersion^ 8 .. 

Several attempts^^! have been made to estimate U for 
Na^CoC^; all of these find that U > W. This is much 
larger than in the organic charge transfer salts such as 
k-(ETJ2Cu 2 (CN) 3 or /3'-[Pd(dmit) 2 ]A , in which U > W 
(Ref. |53 and references therein). 

It is worth keeping in mind the limitations of this 



Hamiltonian as a model for Na^CoO^- For instance, 
the charge ordered state observed in a narrow region 
around x — 0.5 seems likely to be due to the response 
of the electronic system to the ordering of the Na atoms 
which is observed even at high temperatures by x-ray 
diffraction. 8 Therefore one does not expect this phase to 
appear in purely electronic models such as those consid- 
ered here. Nevertheless, before attempting to understand 
more complicated models with realistic band structures 
(including possibly multiple bands) or including the ef- 
fects of Na ordering, it is important to have a firm under- 
standing of strong correlations in simple frustrated mod- 
els. Therefore, despite the remarks above, we present 
below a DMFT study of the Hubbard model on a trian- 
gular lattice Q). In section IVIII Bl we discuss extensions 
to the Hubbard model which allow one to investigate such 
effects. 



C. Kinetic energy frustration 

Geometrically frustrated antiferromagnetism has been 
studied extensively and the triangular lattice provides a 
model system to investigate. Three of the most widely 
used quantitative measures of frustration in antiferro- 
magnets are (i) the degeneracy of the ground state, (ii) 
the magnitude of the entropy at low temperatures, and 
(iii) the ratio of the ground state energy to the total en- 
ergy of maximising the individual interaction energies i^i 

In a non-interacting electron model the only proposal 
we are aware of for a quantitative measure of the geomet- 
rical frustration of the kinetic energy is due to Barford 
and Kim. — They noted that, for t > 0, an electron at 
the bottom of the band does not gain the full lattice ki- 
netic energy, while a hole at the top of the band does. 
They suggested that for tight binding models a measure 
of the frustration is then A = \e™ ax \ — \ e™ m \, where e™ ax 
and e™ m are the energies (relative to the energy of the 
system with no electrons) of the top and bottom of the 
band respectively. This frustration increases the density 
of states for positive energies for t > (negative energies 
for t < 0) which represents an increased degeneracy and 
enhances the many-body effects when the Fermi energy 
is in this regime. However, perhaps a simpler measure 
of the kinetic energy frustration is W/2z\t\, where W is 
the bandwidth and z is the coordination number of the 
lattice. The smaller this ratio, the stronger the frustra- 
tion is, while for an unfrustrated lattice W/2z\t\ = 1. 
For example, on the triangular lattice the kinetic energy 
frustration leads to a bandwidth, W — 9\t\, instead of 
12|t| as one might naively predict from W = 2z\t\ since 
z = 6. 

We show below that geometrical frustration of the ki- 
netic energy is a key concept required to understand 
the properties of the Hubbard model on the triangular 
lattice. In particular it leads to particle-hole asymme- 
try which enhances many-body effects for electron (hole) 
doped t > (t < 0) lattices. 
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It should be noted that geometrical frustration of the 
kinetic energy is a strictly quantum mechanical effect 
arising from quantum interference. This interference 
arises from hopping around triangular clusters which 
will have an amplitude proportional to t 3 which clearly 
changes sign when t changes sign. In contrast on the, 
unfrustrated, square lattice the smallest possible cluster 
is the square the amplitude for hopping around a square 
is independent of the sign of t as it is proportional to 
t A . Barford and Kim£i noted that the phase collected 
by hopping around a frustrated cluster may be exactly 
cancelled by the Aharonov-Bohm phase associated with 
hopping around the cluster for a particular choice of ap- 
plied magnetic field. Thus a magnetic field may be used 
to lift the effects of kinetic energy frustration. The quan- 
tum mechanical nature of kinetic energy frustration is in 
distinct contrast to geometrical frustration in antiferro- 
magnets which can occur for purely classical spins. 



D. The Hubbard model on a triangular cluster 

1. Non-interacting case (U = 0) 

In order to gain some basic insights into some of the 
rich physics associated with the triangular lattice Hub- 
bard model, we now review the exact solution on a single 
triangle. This is the simplest possible model with ki- 
netic frustration and strongly correlated electrons. This 
model already contains some of the features found in the 
solution of the Hubbard model on the triangular lattice 
in the thermodynamic limit and gives significant insight 
into that problem. In particular the solution of this toy 
problem illustrates why the sign of t is important and 
why ferromagnetism arises in the solution to a problem 
that only contains kinetic, local and (via superexchange) 
antiferromagnetic interactions. 

Let us begin by recalling the solution of the two site 
cluster with U — and spinless fermions. Let us label 
the Wannier functions on the two sites |1) and |2). For 
t > the single electron ground state is the bonding 
orbital with wavefunction = + |2)) and energy 

E = —t = —\t\. For t < the single electron ground 
state is the antibonding orbital with |\&) = ^Jg(|l) — 1 2)) 
and E = t = —\t\. (Note that in general changing the 
sign of t reverses the ordering of the energy levels.) Thus 
W = 2t for the two site cluster which, as z = 1 for the 
two site lattice, fits with the general expectation that, for 
an unfrustrated lattice, the bandwidth is given by 2zt. 

The ground state of the U — triangular cluster has 
some significantly different properties compared to the 
two-site cluster. For t > the single electron ground 
state is completely bonding with \^f) = -^=(|1) + |2) + 

|3» and E = —2t = -2\t\, but for t < we cannot 
construct a completely antibonding solution and so the 
single electron ground state is degenerate with |\&) = 
1) - 2|2) + |3)) or |*) = i (|1) + |2) - 2|3)) and 
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FIG. 2: Energy cost in spin polarizing the three site triangu- 
lar cluster with four non-interacting electrons. For t > the 
ground state is degenerate and there is no energy cost in flip- 
ping one spin. In contrast, it costs 3|t| to polarize the t < 
cluster. This indicates that the tendency to ferromagnetism 
of the electron doped t > triangular cluster is much stronger 
than that of the cluster with t < even at the one-electron 
level. 



E = t = —\t\. Thus even at the non-interacting level 
the bandwidth is reduced (W — 3t, z = 2 which implies 
W/2zt = 3/4) by the effects of geometric frustration. 

Moving to noninteracting spin 1/2 fermions, one can 
already see how the three site cluster with t > is easier 
to magnetize than the t < case. In Fig. [3] the energy 
level structure of the t < and t > three site clusters 
containing four electrons are displayed. It costs no energy 
to flip one spin in the cluster, AE = -B(TTTI)--E ; (TITI) = 
0, for t > whereas there is an energy cost of AE = 3\t\ 
associated with flipping a spin in the t < case. We will 
see in section I VII that this tendency to ferromagnetism 
persists in the non-interacting triangular lattice. 



2. Interacting case (U > 0) 

For a strongly interacting system the filling factor n 
will play an important role in the solution. A non-trivial 
case in the triangle relevant to Na^CoC^ is n — |. The 
n = | case is trivial as is simply the case of a single 
(and therefore non-interacting) hole, while the half filled 
case will not be studied in detail in this paper. The 
relevant basis set is rather different from that for non- 
interacting particles and we therefore illustrate the solu- 
tion in Fig. To represent this basis set more precisely 
we will adopt the notation \a, (3, 7) where the values of 
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a, (3 and 7 indicate the state on the first, second and 
third sites respectively. We allow a, (3 and 7 to take the 
values 0, indicating an unoccupied site, f, indicating a 
site singly occupied with a spin up electron, J., indicating 
a site singly occupied with a spin down electron, and 1 [ 
which indicates a doubly occupied sitei^ For t < we 
find that for any finite U the ground state has energy^ 



U>0, t>0: 



(3) 



and is a superposition of singlets 



I*-) = -sinf 
A 



+ — cos( 
A 



I U,U,o) + | TI,o,U) + |o,T|,TI) 
(l 44,4,4}- 1 44,4,4} 



+ (l T,TW>-U>U,T>) 

+(i ra,u)-u,r,u)) 



(4) 



where 



tan6> 



U + Qt- V36t 2 - AtU + u 2 



U-6t+ V36< 2 - AtU + U 2 



(5) 



and the normalisation factor is A = v3 + 3 cos 2 9. 
Clearly, 9 — > as U — > 00 and 9 — » — 7r/4 as {/ — > 0. 
We sketch this wavefunction in the lower part of figure 01 
is a superposition of the resonating valence bond 
state and the states with two doubly occupied sites. The 
fact that 9 — ► as U -^00 indicates that the ampli- 
tude for having two doubly occupied sites vanishes in this 
limit as one expects. Clearly \^ -) has no net magneti- 
zation, this state leads to the paramagnetic state found 
in the thermodynamic limit for t < (c.f. figure [T]l . 
I^-) has significant short-range antiferromagnetic spin 
fluctuations, which are not captured in the purely local 
DMFT treatment that follows, although it is likely that 
these persist in the true thermodynamic ground state (see 
section III A|) . 

For any finite U and t > the ground state has energy 
E+ = —2t + U and is threefold degenerate. All three 
states are spin 1 and three degenerate states correspond 
to S z = —1,0 and 1. The respective eigenstates corre- 
spond to superpositions of triplets and are 



^=1 TI,I,I)-|I,U,I) + U,I,TI), (6) 
1 

7! 



U,U) + I U,U) 



t,TU) + l i,U,T> 



+ (l T,I,TI) + U,T,TI) 



(7) 



and 



V3 



I TI,T,T>-| T,TI,T) + I T,T,U>. (8) 






U>0, t<0: 



cos(9)i 
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FIG. 3: (Color online). A pictorial representation of the 
ground states of the Hubbard model on a triangular clus- 
ter with four electrons. Doubly occupied sites are indicated 
by the presence of two arrows, one pointing up and the other 
pointing down, while unoccupied sites have no arrows. In the 
cases where we have only one doubly occupied site the re- 
maining sites with only one electron can form either a singlet 
(indicated by a thick red line with arrows on the two sites 
pointing in opposite directions) or a triplet (indicated by a 
thick blue line with arrows on the two sites pointing in the 
same direction). Thus, the middle left triangle indicates the 
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state -s(IU,t)i) — MLiit))- The top left triangle indicates 



one the three states | 1,1,11), ^(1 4,1,10 + 1 4,4,U» or 
4, 4, 11}- To keep the figure as simple as possible we have 
not included normalization factors in the figure, but the cor- 
rect normalizations are given in Eqs. HJ, ®, Q an d ©- 
The lower part of the figure illustrates |^-), the ground state 
for any non-zero value of U and t < [c.f. eqn. JIJ] (the 
ground state is degenerate for t < and U = 0). The upper 
part of the figure illustrates all three degenerate states ), 
j* ( 4) and |*|) as these only differ by the value of the S z = 
projection of the triplet. These states are the ground states 
for t < and any value of U [c.f., eqns. JBJ, J7J and @]. |$_) 
depends on the variable which is a function of the ratio U /t, 
given by eqn. ©. 8 — * as U — » oo; 8 — n/4 for U = 0. 



These wavefunctions are illustrated in the upper part of 
figure |2| 

|^+), an d |^+) display short range ferromagnetic 
fluctuations for all finite U, in particular the total spin of 
the system is a maximum along one axis as the ground 
state consists of a superposition of triplets. As the ferro- 
magnetism is associated with kinetic energy frustration, 
we have a Nagaoka type ferromagnet for any finite value 
of U. In the calculations for the infinite triangular lattice 
with t > presented below we will again see this type of 
magnetism for some fillings, but with a phase transition 
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from a Curie- Weiss metal to a Nagaoka ferromagnet at 
some finite value of U, as is shown in Fig. 2] 

In order to quantify the nature and magnitude of the 
spin interactions induced by the correlations we evaluate 
the expectation value of Si- Sj for neigbouring sites on the 
triangle. We find that for t > 0, this expectation value is 
1/12, whereas for t < it is - cos6' 2 /[3(l + cos6' 2 )] which 
becomes more negative with increasing U and tends to 
-1/6 for U 3> \t\. For reference, for the two site Hub- 
bard model, with two electrons, this expectation value 
is always negative, regardless of the sign of t, indicating 
antiferromagnetic spin correlations. Here, we see that for 
the triangular cluster, the sign of t determines whether 
the nearest neighbour spin correlations are ferromagnetic 
or antiferromagnetic. 



III. DYNAMICAL MEAN-FIELD THEORY 

A. Relevance of DMFT to the two dimensional 
triangular lattice 

Before presenting our results an important issue to ad- 
dress is the relevance of DMFT to model {Q which de- 
scribes a two-dimensional system. DMFT is only exact 
in infinite dimensions or for a lattice with infinite coor- 
dination number. We will see below that an important 
role of the triangular lattice, as compared to other, non- 
frustrated, lattices, is not just that the shape of the bare 
DOS changes (which is extremely important) but also to 
make DMFT a good approximation. 

DMFT has become an important tool in the descrip- 
tion of strongly correlated systems. It has provided a 
realistic description of transport and dynamical prop- 
erties of materials such as transition metal oxides and 
layered organic superconductors with frustrated lattices 
such as k-(BEDT-TTF)2X where X is an anion, e.g., 
I 3 or Cu 2 (CN) 3 . The k-(BEDT-TTF) 2 A family have 
a phase diagram (as a function of pressure, uniaxial 
stress, or chemical substitution) in which a supercon- 
ducting phase is in close proximity to a Mott insulating 
phase. In the metallic phase of the k-(BEDT-TTF) 2 A 
there is a temperature scale, T*, at which there is a 
smooth crossover from a Fermi liquid to an incoherent 
"bad metal" [characterized by the absence of a Drude 
peak in the frequency dependent conductivity, a resis- 
tivity of the order of the Ioffe-Regel-Mott limit (H/e 2 a, 
where a is the lattice constant), a thermopower of the 
order of /ce/|e|, and a non-monotonic temperature de- 
pendence of the thermopower and resistivity] The 
temperature dependence of transport properties displays 
a crossover from coherent Fermi liquid behavior to inco- 
herent excitations. In the case of the organic supercon- 
ductors it has been found that the DMFT of the Hubbard 
model gives both a qualitative^ 3 , and quantitative^ de- 
scription of this crossover from a Fermi liquid to a "bad 
metal" in which there are no quasi-particles. 

The crossover temperature scale is related to the de- 



struction of Kondo screening and Fermi liquid behavior 
with increasing temperature (above the Kondo temper- 
ature Tk) in the Anderson model (which is the effec- 
tive model solved in the DMFT equations as described 
above). In the Anderson model the conduction electrons 
are strongly scattered by a (localised) magnetic impu- 
rity for T > Tk- But for T < Tk a singlet forms be- 
tween the impurity and the conduction electrons. In the 
DMFT of the Hubbard model for T > T* the electrons 
are quasi-localised and the electrons on the single site 
treated exactly strongly scatter those in the bath. How- 
ever, for T < T* transport is coherent and the electrons 
only scatter one another weakly, thus Fermi liquid be- 
haviour is regained. This is why, for example, the tem- 
perature dependence of the resistivities of the Anderson 
and Hubbard models are so similar. 

The success of DMFT in describing the transport prop- 
erties and the phase diagram of many organic charge 
transfer salts down to temperatures of about 10-20 K 
(where, for example, superconductivity becomes impor- 
tant) has been rather puzzling given that these materials 
are quasi-two-dimensional and DMFT is only expected 
to be a good approximation in the limit of high dimen- 
sion or co-ordination number. However, the applicability 
of DMFT to low-dimensional systems with large frustra- 
tion is consistent with the fact that for frustrated mag- 
netic models a Curie- Weiss law holds down to a much 
lower temperature than for unfrustrated models 60 - 61 ' 62 
indicating the presence of well formed local moments. 
Recently, Zheng et alm2> calculated the temperature de- 
pendence of the magnetic susceptibility of the antiferro- 
magnetic Heisenberg model on an anisotropic triangular 
lattice. They found that for models close to the isotropic 
triangular lattice (i.e., models with significant magnetic 
frustration) that the Curie- Weiss law (which is a mean- 
field, single site approximation) holds down to relatively 
low temperatures. Deviations from Curie- Weiss behavior 
result from spatially dependent correlations. Hence, we 
expect that a DMFT treatment of the Hubbard model on 
the triangular lattice will be a good approximation down 
to much lower temperatures than unfrustrated models. 
Furthermore, in the "bad metal" region magnetic prop- 
erties such as the uniform susceptibility and spin relax- 
ation rate, can be described by the Heisenberg model 
because the electrons are essentially localized due to the 
proximity to the Mott insulating phase. This means that 
the susceptibility can be fit to a Curie- Weiss form down 
to temperatures much less than the exchange energy J. 
Furthermore, the spin correlation length of the antiferro- 
magnetic Heisenberg model increases with temperature 
T much more slowly for the triangular lattice than the 
square latticed Specifically, at T = 0.3 J the spin corre- 
lation length is only one lattice constant for the triangu- 
lar lattice. In contrast, for the square lattice the correla- 
tion length is about 50 lattice constants, at T = 0.3 J. 63 

The above arguments have been recently tested by 
means of cluster DMFT calculations. These calculations 
show how for the isotropic triangular lattice the solution 



8 



coincides with single site DMFT (in particular a quasi- 
particle peak appears at the Fermi energy). However, as 
soon as frustration is released a pseudogap opens up in 
the one-electron spectra as a result of short range anti- 
ferromagnetic correlations^ 

A further hint that DMFT is a better approximation 
on the triangular lattice than on the square lattice comes 
from the fact that we find that at half filling our calcula- 
tions predict that the Mott transition occurs at U c « 15|t| 
(see Figs. nand|§J). This compares with exact diagonal- 
ization studies on 12 site lattices^ which find that the 
Mott transition takes place at U ~ 12|t|. On the square 
lattice it is known that perfect nesting means that the 
ground state is a Mott insulator for any finite U. How- 
ever, DMFT predicts^ that U c « 12|t| unless antifer- 
romagnetism is included. Thus (without including an- 
tiferromagnetism) DMFT gives a qualitatively incorrect 
result for the (unfrustratcd) square lattice, but a qualita- 
tively correct result for the (frustrated) triangular lattice. 



1. Magnetic susceptibilities 

The magnetization of the system under a small mag- 
netic field, h, is computed from 



m = Mb(«t - n l)> 



(12) 



where raj and raj are obtained from the DMFT solution 
at self-consistency. 

The uniform (q = 0) susceptibility is then obtained 
numerically 5 - 8 from 



X(T) = lim 

h->o oh 



(13) 



It is interesting to compare the uniform susceptibility 
with the local susceptibility obtained from 



=CT) = $>(<!) = A f\s z (0)S z (r))dr (14) 
Jo 



B. Formalism 

The electronic and magnetic properties of the Hub- 
bard model on the triangular lattice Q are analyzed 
below using DMFT^i Previously, various properties of 
doped Mott insulators on an hypercubic lattice have been 
explored within DMFT in the context of the high-T c 
superconductors^ using Quantum Monte Carlo (QMC) 
techniques. Here, we apply DMFT to a frustrated lattice 
and use exact diagonalization and Lanczos techniques to 
solve the associated Anderson impurity problem^ at fi- 
nite and zero temperature, respectively. 

We will now briefly describe the DMFT formalism fo- 
cusing on the relevant equations and quantities of in- 
terest in this work. For a more detailed discussion see, 
for example, Ref. DMFT treats the quantum dy- 

namics on a single site exactly and the remaining lattice 
sites provide a bath with which this site interacts. One 
may map the Hubbard model onto an Anderson single- 
impurity model which must be solved self-consistently. 
The bath is described through the Weiss field, Go CT (*^n)- 
The iterative procedure starts by solving the Anderson 
model for a given choice of Go a (iuj n ). From the on-site 
Greens function, G a (iu> n ) we may obtain the self-energy 
of the system 

(9) 

(10) 
(11) 



which is used to describe the lattice propagator 
G a (iui n ) = y^G^k, ioj n ) 



E 



i 



where /i is the chemical potential, ui n — (2n + 1)tt//3 is 
a Matsubara fermionic frequency and (3 = 1/kgT. The 
above procedure is repeated until a self-consistent solu- 
tion for the lattice propagator is found. 



where S z (t) — raj(r) — raj(r). 

The local susceptibility in the imaginary frequency axis 
is given by 

XlociwJ = I*b e-^ T (S z (0)S z (r))dr. (15) 
Jo 

The frequency-dependent local magnetic susceptibility 
is related2i to the nuclear magnetic resonance (NMR) 
Knight shift K(T) by 



K(T) = A lim Rexzoc(w + if]) 

UJ— »0 



(16) 



when the hyperfine interaction A form factor is indepen- 
dent of wavevector, i.e. in the local limit appropriate to 
DMFT. The NMR relaxation rate l/T x is given by 



1 1 Imxioc(u + if]) 
oc nm 

T\T w^O 7T UJ 



(17) 



where w is a real frequency and rj is an arbitrarily small 
real number. 



2. Numerical methods 

Due to the rapid convergence of the solution with the 
bath siza 5 ^ it is sufficient to use a discrete set of states of 
n s — 6 — 10 sites to model the electronic bath in order to 
obtain reliable results. In what follows we use Lanczos 
diagonalization on n s = 8 sites to calculate zero tem- 
perature properties. In Sec. IIVI we also present results 
for high temperatures for which the Lanczos technique is 
not adequate, we therefore use exact diagonalization to 
solve the Anderson impurity problem. However, as exact 
diagonalization is more computationally expensive than 
Lanczos diagonalization we are limited to n s =6 for these 
calculations. Further it is well known^i that for exact di- 
agonalization there is a low energy scale T ns = E\ — Eq, 
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where Eq and E\ are respectively the ground state en- 
ergy and the energy of the first excited state of the sys- 
tem (recall that there is a finite energy gap to the first 
excited state because of the finite bath size) , below which 
the results of exact diagonalization calculations are not 
trustworthy. We therefore present exact diagonalization 
results for finite temperatures with T > T ns and Lanc- 
zos data valid in the T — > limit. Combining the two 
methods allows us to have a rather complete description 
of the T-dependence of various properties. 



IV. TEMPERATURE DEPENDENCE OF 
MAGNETIC SUSCEPTIBILITY 

We now analyze the behavior of the uniform and local 
magnetic susceptibilities obtained from DMFT for the 
Hubbard model on the triangular lattice for both t > 
and t < 0. We find that the two different signs of the 
hopping integral lead to very different magnetic responses 
(c.f. Fig. 2]) and that Curie- Weiss metallic behavior ap- 
pears when t > but not for t < 0. 




ItlT 



FIG. 4: (Color online.) The imaginary time local spin 
autocorrelation function for both t > and t < with 
n = 1.3, U = 10]*] and (3 = 3/]*]. The larger values of 
(Sz(Q)Sz(t = p/2)) for t > than for t < indicates the 
greater degree of localization of the electrons for t > 0. 



A. Local Magnetic susceptibility 

The degree of localization of the electrons in the trian- 
gular lattice with different signs of t can be explored by 
computing the imaginary time, r, local spin autocorrela- 
tion function, 



H(t) = (S z (0)S z (t)). 



(18) 



Such a correlation function can be related to a spectral 
density A(w) byS 



1 f°° 

H(t) = — / dwexp(-wT)- 



If A{lo) = -A{-w) then 



H(r) 



1 

2^ 



A(w) [e 



-(P-r)u 



(19) 



(20) 



It then follows that H(t) — H((3 — r), as one expects for 
the correlation function of commuting operators^ and 
H(t) should be symmetric about r = (3/2. The function 
H(t) is related to the local susceptibility at Matsubara 
frequencies by equation (TSJ. 

The simplest possible form that the frequency depen- 
dent local susceptibility x; oc (w) [compare equations i|16|) 
and (|17|l ] can take is a "Drude" type form with a sin- 
gle relaxation rate, T, xioc 

(uj) = xo/(l + ioj/T), leading 
to a Lorentzian form for the spectral density. We note 
that for a single impurity Anderson model it was found 
that this Drude form of the relaxation rate is a good 
approximation for temperatures larger than the Kondo 
temperature**! 



In a Fermi liquid, H(t) ~ 1/r 2 for 1/r ~ 1//3 < 
T*. In a Mott insulator, H(t) decays exponentially to a 
non-zero value for long imaginary times r ~ /3/2^i The 
"bad metal" regime represents intermediate behaviour, 
and the value H{t = (3/2) is a measure of the extent to 
which the electrons are localised. 

In the limit in which (S z (0)S z (t)} becomes constant 
one obtains from Eq. (|14fl 



Xioc{T) 



Mb 
T ' 



(21) 



recovering the Curie law for localised magnetic moments. 
In Fig. Hwe have plotted (S z (0)S z (t)) for U = 10|t|, n = 
1.3 at inverse temperature of (3 = 3/|t|, for both t > 
and t < 0. Indeed we find that electrons for the t > 
triangular lattice (weak T-dependence) are more localized 
than those in the t < case (strong T-dependence). Both 
curves are symmetric about r = (3/2 as is required. 

We find that (S z (0)S z (t)) is not constant for either 
sign of t which reflects spin relaxation processes and the 
fact that electrons are not completely localized. The ac- 
tual behavior at large temperatures is not the pure Curie 
behavior of Eq. (|21() but rather, to a good approxima- 
tion, can be fitted to the Curie- Weiss form 



Xioc(T) 



T + T*(n) : 



(22) 



where, T*(n), is again the coherence scale. 

In Fig. Elwe plot l/x; oc (T) from the actual numerical 
results obtained from DMFT for different electron occu- 
pation comparing t > with the t < case. Over a 
broad temperature range the temperature dependence is 
consistent with a Curie- Weiss form. However, the coher- 
ence scale is much larger for t < than for t > 0. 
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FIG. 5: (Color online.) The temperature dependence of the 
inverse of the local spin susceptibility for different electron 
occupations, n. We present results for U = W\t\ and for both 
t < and t > 0. The lines are guides to the eye. For the 
t < cases (only n = 1.3 is shown for clarity) the coherence 
temperature T* <C \t\, whereas for the t > cases T* ~ \t\. 



B. Uniform Magnetic susceptibility 

In Fig. we show the temperature dependence of 
x(T)/%o for the triangular lattice with t > and U — 
10 1 £ | , where xo is the non-interacting uniform susceptibil- 
ity of the triangular lattice. We find that the temperature 
dependence of the uniform magnetic susceptibilities are 
very different for the different signs of the hopping in- 
tegral. The magnetic susceptibility displays Curie- Weiss 
behavior at large temperatures when t > whereas for 
t < 0, x{T) displays the Pauli paramagnetism character- 
istic of a weakly correlated metal. 

We now discuss the t > case in more detail. To be 
specific, for t > and at high temperatures, the uniform 
susceptibility can be fitted to the expression 



X(T) 



T + 6{n) 



(23) 



where 9(n) and /i depend on U and n. Note that the 
temperature scale 6{n) reduces to 9{n) — —Tc, where 
Tc is the critical temperature at which the ferromagnetic 
transition occurs for certain values of n, (1.35 < n < 1.6 
for U — 10|t|; see section fvTI for details). The Curie- 
Weiss behavior described by Eq. (|23[) is valid for large 
temperatures, T > T*(n). This is because in this tem- 
perature regime the electronic system behaves like a set 
of quasi-localized magnetic moments. For T < T*(n) 
the Curie- Weiss behavior crosses over to a behavior more 
reminiscent of a renormalized Fermi liquid with a weak 
temperature dependence. Fitting our results to Eq. 
(EU for T > T*{n), gives 0(n = 1.7) w -0.17|*| and 
6{n = 1.3) w 0.01 \t\ for n = 1.3. For the ferromagnetic 
metal appearing for, say, n — 1.38, the transition to a fer- 
romagnetic metal occurs at Tc = —9(n = 1.4) w 0.05|t|. 



FIG. 6: (Color online.) Curie- Weiss versus paramagnetic be- 
havior of the uniform susceptibility, x(T), for the Hubbard 
model on the triangular lattice with U = 10|t|. The in- 
set shows x(T)/xo (where xo is the zero temperature non- 
interacting magnetic susceptibility) comparing the strong In- 
dependence for t > with the nearly temperature indepen- 
dent behavior for t < (Pauli paramagnetism). For n = 1.4 
ferromagnetism appears below a critical temperature, Tc- 
The main panel displays the inverse of the uniform suscep- 
tibility, xo/x(T), showing the Curie- Weiss T-dependence for 
t > 0. The arrow indicates the coherence temperature T* 
found from fitting the local susceptibility data in Fig. E|to 
Eq. SZZl . 



The effective magnetic moment, \i obtained from the fit- 
ting varies between /i = 0.7/is for n = 1.3 and fi = OAfiB 
for n = 1.7. 

As n is increased one expects that the average moment 
will be roughly proportional to the density of unpaired 
electrons, 2— n, at each lattice site. Thus p? is suppressed 
with increasing n and fixed U . Hence one expects that 
/x 2 varies between fj, — /j,b for n — 1 (half-filled system) 
and /j, = for n = 2 (filled band). This is consistent with 
the values of /x obtained from the fit of Eq. I|23|l to our 
numerical results (see the inset of Fig. 0). As the temper- 
ature is decreased the susceptibility changes its behavior 
so that Xo/x(^) becomes less dependent on temperature 
(this is most obvious for n = 1.7). A crossover from 
Curie- Weiss behavior at high temperatures to Fermi liq- 
uid behavior at low temperatures is encountered. The 
properties of the Fermi liquid state (e.g., effective mass 
renormalization, susceptibility enhancement) at low tem- 
perature will be discussed in the following section. 

In Fig. [7] we show a plot of both temperature scales 
obtained from fitting our data to the Curie- Weiss laws 
of Eqs. ffify and to obtain, T*{n) and 6(n), respec- 
tively. T*(n) increases with the electron occupancy. This 
is to be expected as the system becomes more weakly 
correlated in the n — > 2 limit, remaining coherent at 
temperatures comparable with the non-interacting Fermi 
temperature. On the other hand, 0(n) can change sign, 
which is related to the effective short range magnetic cou- 
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interacting metal is given by - 



FIG. 7: (Color online.) Electron occupation dependence of 
the coherence scale, T*(n), the Curie- Weiss scale 9(n) and the 
effective moment in the uniform magnetic susceptibility, p for 
fixed U = 10|t| and t > 0. These parameters are extracted 
by fitting the results of our DMFT calculations of magnetic 
susceptibilities to Eqs. 1221 and 1231 . For these values of U 
and t ferromagnetism is observed in the range 1.3 < n < 1.6 
(c.f. Fig. 0. 0(n) > for n < 1.3 and no ferromagnetism is 
found at these dopings. p decreases to zero as n — > 2 as the 
density of unpaired electrons decreases. This appears to be 
responsible for the absence of ferromagnetism for n > 1.6. 



pling, J(n), present in the uniform susceptibility, x(T'), 
that does not appear in xioc(T). This can be better un- 
derstood by splitting 6(n) intoSi 



6{n) ~ T*{n) + J(n). 



(24) 



As the system is driven closer to the Mott insulating 
state: n — > 1, 9(n) — > J(n) as T*(n) — > which indicates 
the metal-insulator transition. In this limit the system 
behaves as a Heisenberg antiferromagnet with antiferro- 
magnetic exchange coupling given by J = At 2 /U, with 
U > \t\. A change of sign in 9(n) happens at about 
n = 1.3 in Fig. \7\ signalling a ferromagnetic exchange 
interaction. Indeed, this is approximately the doping at 
which the ferromagnetic state is found for U = 10|t| (c.f. 
Fig. P). Thus, as the electronic occupation increases, 
8(n) changes from positive (antiferromagnetic) to nega- 
tive (ferromagnetic) when J < — T*. This is the thresh- 
old ferromagnetic interaction for the "bad metal" to be- 
come ferromagnetic. At larger doping values the absolute 
yields of the magnetic moment, p, become suppressed as 
the number of unpaired electrons is reduced becoming 
p — > as n — > 2 which coincides with the destruction of 
ferromagnetism. 

In contrast to the rather unconventional metallic state 
found for t > 0, the magnetic susceptibility for t < is 
Pauli-like, moderately enhanced by many-body effects. 
This can be observed in the inset of Fig. |B| where a weak 
temperature dependence for all n (for clarity only n = 
1.3 is shown) is found. To be specific, for T <C £f the 
low temperature behavior of the susceptibility for a non- 



X(T) 
Xo 



p"{e F ) p'{e F f \ 7T- 
p(e F ) p{e F f) 6 



■k 2 B T 2 +0{T 4 ), (25) 



where p(e) is the density of states per spin and the 
primes indicate derivatives of the DOS with respect to e. 
The non-interacting susceptibility at zero temperature is 
given by xo = 2/i^p(e_F). For KOa direct evaluation of 
the term proportional to T 2 gives nearly zero for electron 
doping the triangular lattice. Our numerical results agree 
with temperature dependence given by Eq. with a 

moderate enhancement of x(T)/xo suggesting moderate 
many-body effects. Hence, we conclude that the t < 
electron doped triangular lattice behaves as a renormal- 
ized paramagnetic metal. 



RENORMALIZATION OF QUASIPARTICLES 
FOR DIFFERENT SIGNS OF t 



DMFT typically predicts a metal-insulator transition 
as U is increased to U/W ~ 1 at half-filling i2i As soon 
as the system is doped away from half-filling, the sys- 
tem becomes metallic^ with electrons having their mass 
renormalized by the Coulomb interaction. We find that 
the two different densities of states corresponding to the 
two different signs of t lead to different renormalizations 
of the quasiparticles although W — 9\t\ is the same for 
both DOS. This result is in contrast to the half-filled and 
unfrustrated cases where the results do not depend on 
the sign of t. This is because for unfrustrated lattices 
particle-hole symmetry is respected, whereas particle- 
hole symmetry is broken on the doped triangular lattice. 
Half-filled frustrated lattices are not particle-hole sym- 
metric, but this asymmetry is suppressed at large U, and 
particularly in the Mott insulating state. 

In order to understand the electronic properties and 
different magnetic behavior obtained for t < and t > 
(above), the imaginary part of the self-energy along the 
imaginary axis, ImE(io;) for t < and t > is computed 
and shown in Fig. [S] for U = 10\t\. The quasiparticle 
weight can be extracted from the slope of the self energy, 



lim 



<9ImE(io;) 



d(iuj) 



(26) 



The increasing slope at low frequencies of ImS(iw) as 
n — > 1 indicates a stronger renormalization of the quasi- 
particles as the system gets closer to the Mott metal- 
insulator transition. This behavior is apparent for both 
signs of t. However, for t > electrons are more strongly 
renormalized than for t < 0. For example, for n = 1.3 
the slope of ImS(io;) as oj — > is steeper by a factor of 
about 2 for t > than it is for t < 0. 
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FIG. 8: (Color online.) Comparison of the imaginary part of 
the self-energy for t > with the same quantities with t < 0. 
Results are reported for U = W\t\ and various electron dop- 
ings. The linear dependence of the self-energy at low energies 
indicates Fermi liquid behavior. The slope near iui — is used 
to extract the quasiparticle weight Z, shown in Fig. |U] 



A. Effective mass 



FIG. 9: (Color online.) Dependence of the quasi-particle 
weight Z on the Coulomb repulsion U/\t\ at various electron 
dopings. The quasiparticles become more renormalized as the 
doping is decreased due to the proximity to the Mott transi- 
tion. For a given value of U, quasiparticles on the triangular 
lattice with t < are less renormalized than for t > 0. In con- 
trast, for the half- filled case (n — 1), we find the same behav- 
ior of Z for both t > and t < with a Mott metal-insulator 
transition taking place at a common value U/W ~ 1.66. 



The quasiparticle weight extracted from the self energy 
as a function of the Coulomb interaction U is shown in 
Fig. |5J Z is compared for t > and t < when n = 1.3. 
The effective mass is defined by m* = nib/Z; where m& is 
the non-interacting band mass. For U = 10|t| wc find for 
t < 0, m* = 1.5mf, whereas for t > 0, m* = 3.4mf,. For 
n = 1 we find m* — 4.0ma independent of the sign of t. 
This is an important point as it shows that for the same 
U/W the two different signs of t lead to different renor- 
malizations of the electrons. Within DMFT this different 
renormalization of the quasiparticles can be understood 
from the fact that for electron doping the DOS close to 
the Fermi energy is larger for t > than for t < (c.f., 
the bare DOS shown in Fig. HJ). 

The different renormalization of the electrons found 
for the different signs of t in the electron doped system is 
in contrast to the very similar renormalization found in 
the half- filled case (n = 1) for which Z decreases rapidly 
and in the same way for both t > and t < 0, lead- 
ing to a very similar critical value of U c « 1.65VF at 
which the Mott transition occurs. This value is simi- 
lar to DMFT estimates obtained for the Bethe lattice2i: 
U c ps (1.5- 1.7) IF. The value obtained U c = 15\t\ should 
be compared with the U c pa 12|i| obtained from exact 
diagonalization calculations for 12 site clusters^ The 
similar critical values appearing at half-filling can be ex- 
plained from the fact that as U — ► W, the spectral den- 
sities at each site are strongly modified by the Coulomb 
interaction washing out the fine details of the bare DOS. 
In contrast, for the doped case, kinetic energy effects are 
important. 

In Fig. ^]]we show the variation of the effective mass 



with n for U = 10\t\. Strikingly the effective mass in- 
creases when the system is doped off half-filling close 
to n = 1 and t > displaying a maximum at about 
n = 1.07. This is in contrast to the t < case, and 
what has been found previously for non-frustrated lat- 
tices such as the hypercubic lattice, for which the effec- 
tive mass decreases as the system is doped away from 
half- filling, reflecting reduced correlations. In fact for 
n > 1.07 the effective masses do decrease with increasing 
n, reaching the non-interacting limit (Z — > 1), as n — ► 2. 
A non-monotonic behavior of the effective mass is then 
found close to half- filling and t > 0. This is in contrast 
with the simple behavior encountered for the hole doped 
(n < 1) Mott insulator with a flat non-interacting DOS. 
In that case, the effective mass behaves as^i 

m* /n%b oc 1/x, (27) 

diverging at the Mott transition. This result can also be 
derived from slave-boson theoryiSi 

The effective mass increase with initial doping is a re- 
sult of the initial increase of the DOS for t > as ej? 
moves closer to the van- Hove singularity. The subsequent 
decrease in the effective mass can be attributed to the 
decrease in the effective number of charge carriers which 
suppresses correlation effects. It would be interesting to 
see whether slave boson theory, with the DOS for the tri- 
angular lattice, could reproduce the non-monotonic de- 
pendence of the quasi-particle weight on doping. (This 
would involve solving equations (8) and (9) in Ref. 68). 
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FIG. 10: (Color online.) The variation of the effective mass, 
m* /mi = 1/Z with doping n for both signs of t for U 
10 / Notice in particular the non-monotonic variation of the 
effective mass with n for t > 0. 



B. Sommerfeld- Wilson ratio 

In order to explore the effect of magnetic exchange in 
the metallic correlated state we analyze the Sommerfeld- 
Wilson ratio, 



Rw = hm 



xCO/xo 

° l(T)/lo 



1 



1 



Fg 



(28) 



where Fq is the Fermi liquid parameter which is a mea- 
sure the proximity of the system to a ferromagnetic insta- 
bility (Ffi = -1 at the instability) and 7/70 = 1/Z. The 
dependence of Rw on the electron occupation factors are 
plotted in Fig. El for both t > and t < 0. The ratio 
Rw is the same at n = 1 for t > and t < but they be- 
have very differently as the occupation is increased. For 
the t > triangular lattice Rw is strongly enhanced sig- 
nalling the proximity to a magnetic instability at about 
n = 1.35 and n = 1.65 (see the phase diagram in Fig. ^ 
for t > 0). For f < 0, Rw displays a maximum close to 
n = 1.5. 
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FIG. 11: (Color online.) The variation of the Sommerfeld- 
Wilson ratio, R w = UmT^o{x{T)/xo)/(l(T)/jo) with dop- 
ing n for both signs of t for U = 10[t| and T = [t|/200. 



In Fig. El a comparison of the spectral density A a (u>) 
with the non- interacting DOS per spin, p(e), is made for 
the case n = 1.5 and U = 10\t\. For these parameters 
ferromagnetism occurs when t > 0. Notice that within 
the ferromagnetic metallic phase the spectral function 
of majority spin subband (|) is found to be renormalizcd 
more strongly than that of the minority spin subband (|). 
The lower Hubbard band for spin up electrons contains 
more spectral weight than the lower Hubbard band for 
spin down electrons. One may think of this in terms of a 
photoemission experiment, clearly there are many more 
spin up than spin down electrons which can be extracted 
from a ferromagnet. Alternatively, one may view this 
effect as originating from the preclusion of the hopping of 
a minority spin hole to a neighboring site by the presence 
of a minority spin hole. 



D. Resistivity and Kadowaki- Woods ratio 

Within the DMFT approximation the conductivity is 
given byj2^ 



C. Spectral density 

The different renormalization of the electrons for dif- 
ferent signs of t is translated onto a different redistribu- 
tion of spectral weight induced by the Coulomb interac- 
tion as the spectral weight suppressed at the Fermi en- 
ergy must be transferred to higher energies. We investi- 
gate this by computing local spectral densities, A a (uj) = 
— ilmG -(oj + ii]), where r\ is an arbitrarily small real 
number and 



dep(e) 



£cr(iw) 



(29) 



kcr 



(30) 



where V is the volume. In the limit of low temperatures, 
T — > 0, the resistivity becomes p w AT 2 , with the A 
coefficient given by 



A = 



h |C S | 

e 2 $(e>)/' 



(31) 



The dimensionless constant / = f°° m — , tS -> sa 

1/12, Cs is the to 1 coefficient of the imaginary part of the 
self-energy, Im E(w) « -|C s |(w 2 + (ttT) 2 ), and $(e» is 
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FIG. 13: (Color online.) The doping dependence of the 
Kadowaki- Woods ratio A/^y 2 for t < and t > for fixed 
U — W\t\. The lattice parameter, a = 2.84A relevant to 
Naa;Co02 has been used. Note the non-monotonic doping de- 
pendence of the ratio and the fact that even at half-filling the 
ratios are different for the different signs of t. 



< 0.2 - 
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FIG. 12: (Color online.) Spectral densities for U = W\t\ and 
n — 1.5 for the triangular lattice with t < and t > 0. The 
non-interacting DOS is also shown (red line) for comparison 
and energies are relative to the Fermi energy. The majority 
(T) and minority (J.) spectral densities for n = 1.5 and t > 
are shown as the system is ferromagnetic. 



given by 

^» = ^E(H) 2 ^>-^ (^) 

k x ' 

In the above equations e ~p = /i — ReS(O) = ep, where ep 
is the Fermi energy of the non-interacting systemSi for a 
given electron occupation n. From the T 2 coefficient of 
the resistivity, A, and the specific heat slope at low tem- 
peratures, 7 = 2Tr 2 k 2 3 A(e~F)/3 — TT 2 kgp(eF)/Z, where 
A(ep) and p(£f) are the interacting and non-interacting 
quasiparticle density of states per spin. We can now ob- 
tain the Kadowaki- Woods ratio, 

A h 9\Cy\Z 2 

7 2 e 2 4$(e F )7r 4 /c|/p(e F ) 2 ' { } 

The doping dependence of the Kadowaki- Woods ra- 
tio is plotted in Fig. ED comparing the t > and t < 
cases. The figure displays the non-monotonic dependence 
of the ratio with electron occupation and also how differ- 
ent band structures can lead to different absolute values 
of the Kadowaki- Woods ratio. 



VI. FERROMAGNETISM 

Nagaoka rigorously proved^ that the Hubbard model 
on a connected lattice^ in the U — > oo limit displays 
ferromagnetism when one hole (electron) is added to the 
half- filled system when t < (t > 0). This is due to the 
fact that in this case, the kinetic energy is minimized if 
all the spins are aligned in the same direction. This rigor- 
ous treatment has not been extended to doping by more 
than one hole and it remains an outstanding problem to 
further understand this interesting phenomenon, 72 

DMFT has proved to be an important tool to describe 
ferromagnetism appearing due to local electronic corre- 
lations. The possibility of Nagaoka or metallic ferromag- 
netism in a hole doped infinite dimensional fee lattice 
has been previously analyzed within DMFT. 7 » More re- 
cently DFT(LDA)+DMFT calculations (where DFT is 
density functional theory and LDA is the local density 
approximation) have provided a realistic description of 
ferromagnetism in Fe and Ni^ Our calculations show 
that for sufficiently large U and t > ferromagnetism 
occurs while the system is still metallic as the spectral 
density at the Fermi energy is always finite for non-zero 
doping (c.f. Fig. I12fl . In Fig. [21 we show the sharp 
ferromagnetic transition in the T — * limit for t > 
obtained from our DMFT calculations. 

Increasing U stabilizes the ferromagnetic region in a 
broader electron occupation range. Ferromagnetism is 
found to be more stable for 1 < n < 1.5 than for 1.5 < 
n < 2, indicating the importance of correlation effects as 
the system is closer to the Mott insulating phase as n — > 1 
and sufficiently large U. The ferromagnetic transition 
is found to be sharp as shown in Fig. El even at the 
lowest U value for n = 1.5 where the van Hove singularity 
occurs in the bare DOS. In contrast for t < where, for 
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FIG. 14: (Color online.) Ferromagnetism in the electron 
doped t > triangular lattice from DMFT. The magnitude 
of the spontaneous magnetic moment is shown as a function 
of the electron occupation n. As U is increased above U c (n) 
the system magnetizes spontaneously stabilizing a ferromag- 
netic metal in a broad doping region. Ferromagnetism does 
not appear in the electron doped t < triangular lattice. 

electron doping, there is no van Hove singularity we do 
not observe ferromagnetism. This is, of course, rather 
reminiscent of a Stoner type instability. 

In order to explore the possibility of having Stoner fer- 
romagnetism we can obtain the critical value for fer- 
romagnetism from the Hartree-Fock (random phase ap- 
proximation, RPA) Stoner condition 

U?(n)p(e F ) = 1, (34) 

where p(cf) is the bare DOS per spin. One would find 
Uc(n = 1.5) = and U^(n — 1) = 5 for t > which are 
much smaller values than the ones obtained from DMFT 
calculations (see Fig. |TJ- Furthermore, from the Stoner 
condition we would also expect a ferromagnetic instabil- 
ity at t/f (n = 1.3) = 5.26 and U^{n = 1.90) = 10 for 
t < 0. Clearly, simple Stoner ferromagnetism is inconsis- 
tent with our DMFT results not only at the quantitative 
but also qualitative level as the Stoner criterion would 
predict ferromagnetism for both signs of t for sufficiently 
large, but finite, values of U. Stoner theory is known to 
overestimate ferromagnetic tendencies^ In the present 
case, it would predict ferromagnetism for t < in con- 
trast to the more sophisticated DMFT treatment. This 
difference can be attributed to the on-site dynamical cor- 
relation effects which strongly redistribute the spectral 
weight of the electrons and that are not taken into ac- 
count in standard mean-field theories. 

Vollhardt et aL 76 have given a nice review of the es- 
sential features of Hubbard models that are favourable 
towards to ferromagnetism. (i) The energy dependence 
of the density of states near the Fermi energy should be 
sufficiently asymmetric, with the DOS being larger to- 
wards the top (bottom) of the band for electron (hole) 
doping. A flat band system with a singular DOS at the 



upper(lower) edge the most favorable for ferromagnetism. 
Then the increase in electronic kinetic energy associated 
with spin polarisation is smaller than for a constant DOS. 
(ii) Strong Coulomb interactions narrow the bands. Less 
narrowing occurs for the polarised case because polarisa- 
tion reduces the effect of correlations. 

Our DOS is not of the flat band type, however, the 
energy cost in completely polarizing the non-interacting 
electron doped system for t > is found to be much 
smaller than for t < 0. Indeed, for an electron doped 
system this is given by 



AE = 2 deep{c)- I deep(e), (35) 

J fj. J fJ.p 



which should be negative if ferromagnetism is stable. The 
chemical potentials, po and p p , correspond to the unpo- 
larized and fully polarized systems, respectively and e max 
is the energy at the upper edge of the band. Inserting 
p(e) for t > in expression l|35|l with n = 1.3 we obtain 
AE = 0.603|i|, to be compared with AE = 1.34|i| for 
t < 0. Interestingly, for the square lattice, AE = 0.8|t|, 
which is between the values of AE for the t < and 
t > triangular lattices. Hence, although unsurprisingly 
ferromagnetism is not stable for the non-interacting sys- 
tem regardless of the sign of t, the ferromagnetic tenden- 
cies are clearly stronger in the electron doped triangular 
lattice for t > than they are for t < in agreement 
with the conclusions derived from the three site cluster 
analysis (sect. Ill Dl and Fig. 0). 

Finally, it is worth comparing previous studies of the 
t — J model on the triangular lattice with our DMFT 
results. A Curie- Weiss metal as well as ferromagnetism 
have been previously found from calculations using high 
temperature series expansions, for the hole doped trian- 
gular t — J models Our results for the electron doped 
lattice are equivalent to the hole doped case discussed 
in Ref. once the sign of t is reversed. More specifi- 
cally, ferromagnetism and Curie- Weiss metallic behavior 
appear in the high temperature expansion for the hole 
doped system with t < while a weak temperature de- 
pendent susceptibility appears when t > 0. Furthermore, 
a heavy fermion state, analogous to our strongly corre- 
lated metallic state for t > 0, is predicted to exist in the 
hole doped system with t < 0. The qualitative agreement 
of the magnetic response obtained from DMFT compared 
with the high temperature studies (which includes both 
local and non-local correlation effects) of the t — J model 
indicate the importance of the local aspects of electronic 
correlations. This has also been pointed out in the con- 
text of metallic magnetism for which a good a description 
of the electronic and magnetic properties of Fe and Ni is 
attained from a local theory such as DMFT. 74 
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VII. COMPARISON WITH EXPERIMENTS ON 

Na^CoOs 

As this work has been largely motivated by experi- 
ments on Na K Co02 we now consider what the lessons 
learned from a DMFT study of the Hubbard model on 
a triangular lattice might have to say about Na x Co02- 
First, our results show that the details of the electronic 
structure, and in particular the DOS, play a crucial role 
in determining the physics of strongly correlated frus- 
trated systems such as Na x Co02- This is something of a 
pyrrhic victory in that this very result tells us that there 
is little hope of quantitative agreement between experi- 
ments on Na K Co02 and calculations based on a simplified 
one band models with nearest neighbour hopping only, 
such as those presented above. Indeed our results do not 
show any such quantitative agreement. 

In particular for the triangular lattice with t < 
(which, of the two cases we considered, gives the band 
structure closest to those suggested by both ARPES and 
LDA-DFT) the DMFT of the electron doped Hubbard 
model closely resembles a weakly interacting metal. This 
is clearly not what is observed in Na^CoC^, even at a 
qualitative level. 

On the other hand many of the results for t > are 
qualitatively consistent with the picture drawn by exper- 
iments on Na K Co02- In particular we propose that the 
experimental 'Curie- Weiss metal' is little more than the 
"bad metal" in the extreme situation of very strong elec- 
tronic correlations and high frustration. Important ex- 
perimental support for this hypothesis comes from the 
fact that a Fermi liquid like resistivity has only been 
observed^ when T<1K which suggests that T* ~ 1 K. 
An obvious objection, that T* predicted by DMFT will 
not be this small for any reasonable parameters, will be 
discussed below (Sec. I VIII A 111 . 

Our results for t > show that magnetism is a gen- 
uinely quantum many-body effect, in particular the mag- 
netism has little resemblance to simple Stoner ferromag- 
netism. We propose that the observed A-type antiferro- 
magnetism results from in plane Nagaoka ferromagnetism 
with a weak interlayer antiferromagnetic coupling. The 
observationiS*i& that the strength of the effective in-plane 
ferromagnetic coupling in Na^Co02 is the same order of 
magnitude as the effective interlayer antiferromagnetic 
interaction despite of the highly two-dimensional crystal- 
and band-structures of Na x Co02 is naturally explained 
in this scenario as the effective ferromagnetic interac- 
tion in the Nagaoka is much weaker than the intrin- 
sic antiferromagnetic interaction which gives rise to the 
phase. This can be seen from Fig. [23 here 6(ri) be- 
comes positive and ferromagnetism is suppressed by re- 
ducing the doping towards the critical doping which is the 
lower bound for which ferromagnetism is observed. The 
intrinsic antiferromagnetic interactions due to superex- 
change are given by J = At 2 /U = 0.4|i| for U = 10|i|, 
which we take in the calculations above. Therefore even 
though J = At 2 /U > Jl = 4t 2 L /f7, where Jj_ and t± 



are respectively the inter-plane exchange constant and 
hopping integral, near the phase transition to ferromag- 
netism J 3> 0(n) and thus the inter-plane antiferromag- 
netic exchange may be the same order of magnitude as 
the in-plane ferromagnetic exchange, as is observed ex- 
perimentally in Naa;Co02. 

To qualitatively test these ideas calculations with a 
more realistic band structure are required. LDA calcula- 
tions predict that the dispersion of the bands is rather dif- 
ferent to the dispersion obtained from a nearest-neighbor 
isotropic triangular lattice. In particular the LDA band 
is much flatter close to the T-point. Indeed this is just the 
sort of change that is likely to have a significant effect on 
the results of the type of calculations we have presented 
above. In the current work we have also neglected the 
possibility of Fermi surface pockets near the Appoints 
arising from the e' 2g band^S^ Recall that, for example, 
we demonstrated that the magnetism and Curie- Weiss 
metal behavior we observed in the triangular lattice de- 
pend crucially on the DOS near to the Fermi level and 
not just on the DOS at the Fermi energy. Therefore, even 
if this band does not actually cross the Fermi energy, as 
some studies have indicated,— it may play an impor- 
tant role in Na a Co02. Interestingly, the e' 2g band can be 
fitted to a t > tight-binding model close to the Fermi 
energy. Hence, the e' 2g bands may play an important role 
in producing the Curie- Weiss metallic state. Anyway, the 
generalization of these results to multiband systems^ is 
likely to be important for a complete understanding on 
the magnetic and transport properties of Naa;Co02. 

To quantify these remarks we now present a compar- 
ison of the DMFT results for the triangular lattice with 
the experimentally measured properties of Naa:Co02. 



1. Electronic heat capacity 

The experiment a&S value of the linear- T coefficient 
of the specific heat (C v « 7T for T < T* with T* = 1 K) 
in Na x Co0 2 and 0.7 < x < 0.82 is 7 = 25 - 30 mJ/(mol 
K 2 ) in zero magnetic field. This corresponds to an ef- 
fective mass enhancement m* /mb ~ 3 — 4 where mt is 
the LDA band mass. On applying an external magnetic 
field of H = 14 T, 7 is suppressed leading to 7 = 0.02 
J/(mol K 2 ) (Ref. UJ§ implying at 20-30% decrease in 
to*. The observed temperature and magnetic field de- 
pendence also has two unusual features. In zero field, 
C V (T)/T is non-monotonic, but becomes monotonic in a 
field of 14 T (see Fig. 3 in Ref. 0). The reduction of 
C V (T)/T as the temperature increases from 1 to 5 K is 
consistent with the destruction of quasi-particlcs in this 
temperature range (compare Fig. 37 in Ref. l2l|) 

Our DMFT calculations for n = 1.7 and U = 10|i|, pre- 
dict a weak renormalization of the electrons, m* /mb ~ 
1.1, for t < and to*/to » 1.5 for t > 0. However, in 
the discussion to follow it will be important to note that 
as the system is driven towards half- filling, n — > 1, we 
find effective mass enhancements which are comparable 
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to or larger than experimental values (see Fig. I10|) . Note 
that effective mass enhancement for t < is rapidly sup- 
pressed (Fig. II C)|t as n is increased above n > 1.2 becom- 
ing effectively a weakly correlated metal for this doping 
region. But, for t > we find that the mass enhancement 
remains large to much higher dopings than for t < (for 
t > the mass enhancement is comparable or greater 
than that observed experimentally for x < 1.5). 



2. Low-temperature resistivity 

Experimentally^ the resistivity varies as, p(T) w po + 
AT 2 , as is expected for Fermi liquid, below T w 1 K. The 
coefficient of the quadratic term is measured^ to be A = 
0.96 pfl cm/K 2 for Nao.7Co02- On applying an external 
magnetic field up to H = 16 T, the quadratic coefficient 
of resistivity is decreased to A — 0.22 pFt, cm/K 2 and T* 
is increased to 4 K. 

DMFT predicts Fermi liquid behavior below T*, the 
low energy coherence scale. The resistivity coefficient 
predicted by DMFT is A = 0.0053 pfl cm/K 2 for t < 0, 
n = 1.7 and U = 10|t|, which is more than two orders 
of magnitude smaller than the experimental value. How- 
ever, at half-filling and U = 10j*| , A = 0.81/ifi cm/K 2 , 
which is comparable to the experimental result. For 
t > 0, n = 1.7 and U = 10\t\ we find A = 
cm/K 2 , which is the same order of magnitude as in exper- 
iments of Na^CdC^. Therefore, either the details of the 
DOS or driving the system towards the Mott transition 
could be responsible for the large value of A observed 
experimentally. However it is important to stress that 
in our calculations the coherence scale, T* s» 100 K for 
U = 10|t| at half-filling, is two orders of magnitude larger 
than the temperature below which p(T) w po + AT 2 in 
the experiments^ We will discuss the reasons for this 
below (see Sec. IVlll Alf) . 



3. Kadowaki- Woods ratio 

The Kadowaki- Woods ratio of Nao^CoOa, is experi- 
mentally founoji 9 - to be A/j 2 = 60ao, with o,q = 10 _5 p:f2 
cm mol 2 K 2 /mJ 2 being the constant value found in the 
heavy-fermion materials. The large values of the ratio 
have been recently discussed by Hussey^S who explained 
that when volumetric rather than molar quantities for 
the specific heat, the values of A/7 2 are comparable to 
the ones in the heavy fermion systems. From the DMFT 
calculations with t < we find that A/7 2 = 64.2ao for 
n = 1.7 whereas at half-filling A/7 2 — 75ao, implying a 
weak dependence with the electron doping. A stronger 
dependence is found from DMFT for t > for which the 
ratios vary between A/7 2 = 100ao for n = 1.7 and 170ao 
at half-filling (see Fig. H3|) . Thus, the absolute values 
depend strongly on the DOS. 



4- Uniform magnetic susceptibility 

Experimentally, the magnetic susceptibility is found to 
have a Curie- Weiss form, ~ l/(T+6), and a magni- 
tude much larger than the Pauli susceptibility expected 
for weakly interacting metal across a large range of dop- 
ings that have metallic ground states&ii Our DMFT cal- 
culations for t < display Pauli paramagnetism, which 
is typical of weakly correlated metals, in qualitative dis- 
agreement with experiments. However, for t > a Curie- 
Weiss metal is observed across a wide range of dopings. 
In the discussion that follows it will be important to recall 
that the DMFT prediction of the susceptibility of a metal 
close to the Mott transition is that it has the Curie- Weiss 
form with 9 = T* + J > 0. Thus, a behaviour similar 
to the experimentally observed Curie- Weiss metal is pre- 
dicted by DMFT for t > [i.e., when p(ep ) is large] and 
near to the Mott transition. 



5. Sommerf eld- Wilson ratio 

The experimentally measured Sommerfeld- 
Wilson ratioi 9 - of Nao.7Co02 is Rw = 
lim T ^ (x( r )/xo)/(7/7o) = 7.8. For the hypercu- 
bic lattice, the Wilson ratio predicted by DMFT goes 
to zero close to the Mott transition, i.e. at half-filling, 
but becomes finite as soon as the system is doped away 
from half-filling. 21 On the triangular lattice with t < 
we find that the Wilson ratio attains its largest values 
at n = 1.6 at which Rw = 1-7 and decreases to 1.3 as 
n — * 2 with U = 10\t\. In the limit n — * 1, the calculated 
Rw = 0.8. On the other hand for t > 0, Rw becomes 
extremely large at dopings close to the ferromagnetic 
instability. For example, for n ~ 1.65, Rw is already 
about 6. One should note that experimentally, x = 0.7 
is extremely close to the critical doping at which the 
A-type antiferromagnetism appears. Therefore the 
large Sommerfeld-Wilson ratio observed is probably a 
direct consequence of the proximity of Nao.7Co02 to the 
magnetic transition. This analysis would also suggest 
that as the doping is decreased from x — 0.7 then Rw 
will also decrease, although calculations with a more 
realistic band structure are required to confirm that 
DMFT makes this prediction. 

6. Effect of a small external magnetic field on the low-T 
resistivity, A, and specific heat, 7, coefficients 

A rapid suppression of the low-T resistivity coefficient 
A and specific heat slope 7 and an increase in T* are ob- 
served when an external small magnetic field is applied 
to Na a Co02. [The details of the experimental observa- 
tions of the behaviour of 7 (Ref. [?7j and A (Ref. Il9t) are 
summarised in sections IVII ll and lVII 2l respectively.) 

For the single impurity Kondo model the temperature 
dependence of the heat capacity has a maximum near the 
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Kondo temperature. In a magnetic field, once the Zee- 
man energy, gpsB, is comparable to the Kondo energy, 
UbTk, this maximum shifts to higher temperatures 
Based on this we would expect that the coherence tem- 
perature becomes larger with increasing field. 

DMFT for the Hubbard model in a magnetic field 
shows that in the metallic phase of a weakly interacting 
metal the quasi-particle weight Z increases with field, 
consistent with the predictions of Stoner theory. (See 
Fig. 5 of Ref. l58j) Consequently, the specific heat coeffi- 
cent 7 will decrease with field. Similarily, the resistivity 
coefficient A will also decrease. In contrast, a strongly 
correlated metal close to the Mott transition displays a 
weak suppression with magnetic field (see inset of Fig. 9 
in Ref. l58i). 



of T* predicted by DMFT are comparable with values of 
9 « 150 K found experimentally^ from the Curie- Weiss 
fit of the spin susceptibility, x{T) ~ 1/(T + 9), at large 
temperatures. 

Spin susceptibility experiments on Na K Co02 show a 
clear departure from the simple high temperature Curie- 
Weiss susceptibility at temperatures of order the coher- 
ence temperature predicted by DMFT (but orders of 
magnitude above the temperature where the resistivity 
is observed to vary quadratically with temperature). Be- 
low this temperature the spin susceptibility is strongly 
enhanced in contrast to the DMFT prediction of a Pauli 
susceptibility below T* . Although on this front one 
should beware of the major differences^ between the sus- 
ceptibilities measured by Foo et al& and those measured 
by Prabhakaran et a/ii 



Failure of one band models with nearest neighbour 
hopping 



Spin fluctuations 



The disagreements highlighted above between the one 
band model with nearest neighbour hopping and exper- 
iments on Naa;Co02, particularly for t < which is the 
appropriate sign of t for Naa;Co02, does not appear to 
be a result of the DFMT approximation. Our results 
are entirely consistent with the series expansions^ and 
RVE 12 results for the t — J model on a triangular lattice 
in this regard. Therefore we believe that simple one band 
models with nearest neighbour hopping only cannot ac- 
count for the observed behaviour of Naa:Co02. This is 
an extremely important result given the number of the- 
oretical papers on both NazCoC^ and Na x Co02 ■ J/H2O 
that are based on this type of model. In the next sec- 
tion we examine what kind of models or approximations 
are required to give an account of the strong correlation 
effects in NazCoC^. 



VIII. MORE REALISTIC MODELS OF Na x Co0 2 
A. Beyond DMFT 

Although DMFT calculations on the triangular lattice 
suggest the proximity of the metal to a Mott insulating 
phase several aspects of the observed low temperature 
behavior of Naa:Co02 are still difficult to understand. 
Foremost amongst these is why is the value of T* esti- 
mated from resistivity so much smaller than the obtained 
from the susceptibility and DMFT? 

The experimental value for the coherence temperature, 
T* = 1 K, as determined from the temperature below 
which p(T) = po + AT 2 , is two orders of magnitude 
smaller than the DMFT result, T* m 0.l\t\ m 100 K 
for U = 10\t\ at half-filling. For t > there is only a 
weak increase in T* with doping (see Fig. |5j). Driving 
the system closer to the Mott transition by increasing U 
would lead to a strong suppression of T* , e.g. T* — > 
as U — > U c at half-filling. On the other hand, the values 



All of the results presented in this paper were ob- 
tained from a purely local theory, DMFT. Therefore 
these results neglect short range ferromagnetic fluctua- 
tions which are likely to become increasingly important 
as one dopes towards the ferromagnetic phase. Indeed 
such ferromagnetic fluctuations have already been ob- 
served by neutron scattering in Nao.7Co02«^*^*^ 

The importance of these fluctuations can be addressed 
within Moriya's 82 spin fluctuation theory of the elec- 
tronic and magnetic properties of nearly ferromagnetic 
metals. The theory assumes that the dynamic spin sus- 
ceptibility has the form 



X (q = O,cj = 0) 



1 



(36) 



in two dimensions. £ is the ferromagnetic correlation 
length and T is a phenomenological parameter. This has 
the same frequency and wavevector dependence as one 
obtains from an RPA treatment of a Fermi liquid metal. 

Within this framework, the electronic properties can 
be sensitive to an applied magnetic field as, at low tem- 
peratures, it suppresses the ferromagnetic fluctuations, 
driving the system into an unrenormalized Fermi liquid 
state. The low temperature scale Tq is defined by 8 ^ 



To 



2vr 



(37) 



where i:q 2 B is a measure of the area of the first Brillouin 
zone. 

The spin susceptibility is enhanced as the temperature 
is lowered departing from the Curie- Weiss law as a result 
of the enhancement of ferromagnetic fluctuations at low 
temperatures. Also a small magnetic field is expected 
to rapidly suppress A and 7 as experimentally observed 
in Na x Co02- The above arguments are consistent with 
evidence of strong ferromagnetic fluctuations acting in 
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the Co-planes for x ~ 0.75 from recent neutron scattering 
experiments^ In order to determine the importance of 
these spatial correlations it is desirable to compare the 
measurements of uS) in Ref. 14 with the form (|30ll . 
If at a frequency u> one observes a peak of width l/£ 
in wavevector space then Tn ~ uj(qB£,) 3 - One sees from 
Figures 1 and 4 in Ref. [Tj that for u> ~ 10 meV that 
qsS, ~ 0.4. Hence, T ~ 1 — 10 K and so the ferromagnetic 
fluctuations could be the origin of the low energy scale 
in Na^OoO^. 



B. Beyond the single band Hubbard model 

1. Role of ordering of the sodium ions 

We believe that the A-type antiferromagnetic phase 
observed experimentally corresponds to the ferromag- 
netic phase predicted by DMFT and many of the metallic 
state properties can be interpreted in terms of their prox- 
imity to a ferromagnetic transition and strong Coulomb 
repulsion effects. However, a number of experiments sug- 
gest that Na^CoC^ behaves like it is more strongly corre- 
lated as the doping, x, is increased. These results are dif- 
ficult to understand from a simple model which only con- 
tains the Co planes as such models become more weakly 
correlated as the electronic doping is increased due to 
the small number of holes left in the Co-planes. Below 
we propose that the ordering of the Na atoms observed 
at many dopings provides a natural explanation of these 
effects as it can effectively drive the system towards half 
filling by introducing a spatial modulation of the site en- 
ergy of the Co atoms. The suggestion that Na order- 
ing plays a crucial role that has not been appreciated 
until now is made more plausible by the very sensitive 
dependence of the magnetic and transport properties on 
the exact details of the DOS that we have demonstrated 
above. 

It seems likely that the charge ordering of the electrons 
in the C0O2 layers observed at x = 0.5 is induced by Na 
ordering*? Recently, experimental studies by Zandbergen 
et alJH and Roger et alJ&, first principles calculations 
by Zhang et alm^ and classical electrostatic calculations 
by Zhang et al. and Roger et al. have shown that as 
x is varied not only is the doping of the C0O2 planes 
varied, but so is the Na ordering. However, there re- 
mains debate over the exact nature of the Na order at 
specific values of xiSi^iS 5 . There is additional evidence 
from NMR experiments of the presence of inequivalcntly 
charged Co ions with Co 3+ (S=0) and Co 4 + (S=l/2) 
induced by the charged Na layer of ionsiZ 9 . This effect 
implies that a large fraction of the electrons in the Co 
planes behave as S = 1/2 localized moments. This is 
precisely the behavior expected close to the Mott insu- 
lating phase in a Hubbard model close to half-filling as 
described by DMFT. 

To make this proposal more concrete we analyse four 
special =0.33, 0.5, 0.71 and 0.75. For each of 



these we take the orderings proposed by Zhang et alM. 
and Zandbergen et alJQ although we note that in two of 
these cases (x — 0.71 and 0.75) Roger et al. have pro- 
posed different Na ordering patterns which involve the 
formation of 'vacancy droplets'. We stress that this work 
has nothing to contribute to the debate over which pat- 
terns are, in fact, realised and we take these patterns at 
x = 0.71 and 0.75 merely to exemplify our point. In the 
cases of x = 0.33 and 0.5 we find that the Na ordering 
introduces two Co sublattices while for x = 0.71 we find 
three Co sublattices and for a; = 0.75 we find four Co sub- 
lattices. For this analysis one needs to consider the two 
possible locations of Na in Na^Co02, referred to as Na(l) 
and Na(2)i2£ The C0O2 planes are at z = 0.25c, where c 
is the lattice constant perpendicular to the cobaltate lay- 
ers. The Na(l) sites lie directly above (at z = 0.5c) and 
below (z = 0) the Co atoms, while the Na(2) sites are 
above or below the centres of the triangles in the triangu- 
lar lattice formed by the Co atoms. The simplest model 
which can include the effects Na ordering is a Hubbard 
model which includes the Co site energy modulated by 
the number of nearest Na sites that are occupied. Clearly, 
one expects that occupied Na(l) and Na(2) will modulate 
the Co site energy by different amounts which we denote 
£1 and £2 respectively. If /i is the chemical potential for 
x = the Hamiltonian is then given by 

<ij>,a i 

where Nu and A^j give the number of occupied nearest 
neighbour Na(l) and Na(2) sites respectively. This model 
assumes that the electronic degrees of freedom relax on 
much shorter time scales than the Na ions move, which 
is reasonable because of the relative masses, even though 
the effective quasiparticles mass is an order of magnitude 
larger than the bare electronic mass. 

To simplify our notation we now introduce the follow- 
ing nomenclature, 





= £2 - 


- /' 




£B 


= 2£ 2 


- M 




sc 


= 3£ 2 


- M 




£D 


= 4£ 2 


- M 




sp 
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- £ 2 - 


■ A* 


e Q 


= ei4 


-2£ 2 


- V 




= EH 


- 3£ 2 


- V 



This allows us to sketch the relevant lattices in Fig. ^| 
It can be seen that for x = 0.33 we have a honeycomb 
lattice of sites with site energy ea with the central sites in 
the lattice having site energy Ep. For x = 0.5 the lattice 
contains alternating chains of sites with site energy eb 
and Ep. With x = 0.71 we find the lattice contains sites 
with site energies Ep> , Ep and eq with all three sublattices 
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x=0.33 




x=0.5 




m y 



V H- 



x=0.75 




FIG. 15: (Color online.) Different sublattices of cobalt ions 
due to commensurate ordering of sodium ions. Pictorial rep- 
resentation of the Hamiltonians given by eqn. 1381 are shown 
for x = 0.33, 0.5, 0.71 and 0.75. The different letters denote 
the different site energies defined by eqn. 1391 . so that, for 
example a P represents a site with site energy ep. These 
different site energies result from the different Na order pat- 
terns, which have been observed experimentally— and pre- 
dicted from first principles calculations!-- The Na ordering 
for the quoted filling factors are taken to be those shown in 
Fig. 2 of Ref. |H Note that Roger et al£± proposed differ- 
ent ordering patterns for x = 0.71 and 0.75, which involved 
'vacancy clustering'. 



forming stripes. For x = 0.75 a complicated arrangement 
of interlaced stripes of sites with site energies £c, £p, eq 
and er is formed. This model l|38l) will clearly make 
very different predictions from the Hubbard model if the 
variations in the site energies caused by Na ordering are 
sufficiently large. For example for x = 0.5 and £p>eg 
(as one might expect) and P~H'a charge ordered in- 
sulator is likely to be the ground state of Hamiltonian 
||SSJ, this is the state observed in Nao.sCoC^. Roger 
et al. have estimated from classical electrostatics that 
Ae = Ep — Eb — 100 meV which is consistent with this 
scenario. Further this kind of effect will drive the sys- 
tem closer to Mott insulating phase if the occupation 
of one or more of the sublattices is strongly suppressed. 
Thus we expect that the DMFT treatment of (SB) will 
have significantly better qualitative agreement with the 
experimentally observed properties of Na^CoC^ than the 
DMFT treatment of the Hubbard model on a triangular 
lattice does. Finally we note that, if one of the sublat- 
tices of x = 0.5 can be 'integrated out' of the effective low 
energy theory due to a large disparity in the site energies 
we are left with a Mott insulator on a rectangular lat- 
tice. Along one side (the horizontal direction in Fig. I15|) 
of the rectangle it can be seen that the antiferromagnetic 
exchange constant is J = At 2 /U . In a similar manner 
superexchange leads to an exchange constant 



J' 
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1 
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U 2Ae + U 2Ae 



(40) 



in the perpendicular direction. Additionally there is a 
(diagonal) next nearest neighbour exchange interaction 
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It is interesting to consider the three limiting cases: 



for U > Ae, 
for U < Ae, 
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fort/ = Ae, y = _ « 1 and - = -. 

In particular note that for U ^> Ae we may have J' ~ J. 
Neutron scattering experiments suggest that the charge 
ordered phase of Nao.sCoC^ has long range Neel order— 
which suggests J ~ J' > J" and is therefore consistent 
with this model when U ^> Ae ^> t, which is precisely 
the regime that the separate estimates of these three pa- 
rameters suggest we are in. 

Before moving on from the effects of Na ordering it is 
interesting to note that this type of problem is actually 
rather general in doped systems. For example, in most 
models of the cuprates it is assumed that the only effect 
of varying the doping is to vary the electronic density 
in the C0O2 layers. However, if it is confirmed that Na 
ordering plays the role we have proposed in Na^OoC^ 



21 



then this may be an important step towards a general 
understanding of the role of dopant impurities in strongly 
correlated systems. 

We stress that the above proposal relies on Na order- 
ing and not disorder. It is well known that disorder only 
plays a significant role in DMFT when there are rather 
extreme levels of disorder^ Therefore one would not ex- 
pect our results to be significantly affected by the low 
levels of disorder typically found in Na^OoC^. However, 
recently the role of disorder on the band structure of 
Na^CoC^ has received some attention^ If disorder does 
induce significant changes in the DOS this would be an 
important factor for the results reported in this paper. 
More generally, the role of disorder in Na^CoC^ • 2/H2O 
will no doubt prove an important probe of whether the 
superconducting state is conventional or not. 88 

2. Multiple bands 

The model, as written in eqn. (|38|l . only contains a sin- 
gle band. As discussed in the above band structure calcu- 
lations suggest a simple one band model is not sufficient 
to describe the band structure of NazCoC^. Although 
the ai 9 band can be roughly fitted to a tight-binding 
model with t < which is consistent with the large hole 
pocket around the T point, additional small hole pockets 
may be present in the Fermi surface associated with the 
e' 2g bands crossing the Fermi energy. These bands can 
be fitted to a t > tight-binding band which is likely 
to favour ferromagnetism. Such a multi-band model is a 
straightforward generalisation of (|38f) giving 

iav 

- Ujvv{c\ V(T c jv , a + H.c.) (42) 

ivv' 

where the labels v and v' refer the different orbitals. A 
multiband DMFT treatment is likely to describe both of 
the following effects: (i) the proximity of the metal to 
a Mott insulating phase leading to renormalized quasi- 
particles and Curie- Weiss behavior for T > T* and (ii) 
the strong ferromagnetic fluctuations appearing at very 
low temperatures. Several models have been recently 
proposed to describe Na^CoO^. In comparison to sin- 
gle bancUSiSii and multiband Hubbard models^^SiS, our 
proposed model explicitly contains the periodic poten- 
tial of the Na layers acting on the Co planes in com- 
bination with the Coulomb interaction. This is an im- 
portant ingredient to understand the unexpected Curie- 
Weiss susceptibility and large effective mass at x — ► 0.7. 
At these large dopings the system is expected to behave 
as a weakly correlated metal. We believe that our model 
is sufficient to fully understand the interplay between 
these two essential ingredients of the problem. 



The problem of applying DMFT to a Hubbard model 
with two inequivalent sublattices has been considered be- 
fore. Chitra and Kotliar 90 considered the Hubbard model 
at half filling in the presence of antiferromagnetic order. 
This leads to a different chemical potential on the two 
sublattices. Bulla and collaborators^ considered charge 
ordering in the extended Hubbard model at one-quarter 
filling. For the case of the Bethe lattice the model was 
mapped to a pair of Anderson impurity models, one for 
each of the two sublattices. For V (the Coulomb re- 
pulsion between electrons on neighbouring lattice sites) 
larger than about t and U — 2t the quasi-particle weight 
(and presumably the coherence temperature also) is or- 
ders of magnitude smaller than the non-interacting value 
and there is a pseudogap in the density of states. Based 
on the above we would expect that for x = 0.5 the ground 
state will be insulator if \ep — Eb\ ~ U > \t\ due to the 
two different sublattices. 



IX. CONCLUSIONS 

We have presented a DMFT analysis of the electronic 
and magnetic properties of the doped isotropic triangu- 
lar Hubbard model. An important result of our work 
is the large effect of the bare DOS on its magnetic and 
electronic properties. In particular, a Curie- Weiss metal 
is found for t > and U > W . In contrast, Pauli para- 
magnetism is found for t < regardless of the magni- 
tude of U . We find a larger degree of localization and 
a larger renormalization of the electrons for t > than 
for t < due to the different DOS at the Fermi energy. 
The spin susceptibility crosses over from Curie- Weiss be- 
havior, at large temperatures, to a renormalized Fermi 
liquid as T — ► 0. At low temperatures, the uniform spin 
susceptibility is strongly enhanced for t > due to the 
proximity to the ferromagnetic transition. This is in con- 
trast to the behavior for t < 0, for which only a small 
enhancement of x(T — > 0)/xo is found. The stronger 
renormalization of the quasiparticles found for t > than 
for t < when the system is doped away from half-filling 
is a consequence of a larger DOS near the Fermi energy 
in the t > case. 

For t > we find a metallic ferromagnetic state when 
U is sufficiently large. The ferromagnetism is signifi- 
cantly different from Stoner ferromagnetism in that it 
has a strong local moment character and the criterion 
for ferromagnetism is not simply related to the DOS at 
the Fermi level, but depends on the features of the DOS 
over a significant energy range around the Fermi level. 
No ferromagnetism has been found for any band filling 
or any value of U when t < 0, which is qualitatively dif- 
ferent from the prediction of Stoner theory. We conclude 
that the ferromagnetism we have observed is much more 
closely related to Nagaoka ferromagnetism. The differ- 
ent behaviour of the two lattices is due to kinetic energy 
frustration which dramatically reduces the energy cost of 
polarizing the system for t > 0. We have therefore pro- 
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posed that the A-type antiferromagnetic phase observed 
in Na^CoC^ fori > 0.75 is the result of in plane Nagaoka 
type ferromagnetism. 

The different behaviors encountered away from half- 
filling for the different signs of t contrasts with the system 
sufficiently close to the Mott insulating phase (n — ► 1 and 
U > U c or U — > U c and n = 1). In this case, the model 
displays very similar properties including Curie- Weiss 
susceptibility with 9 > (due to antiferromagnetic su- 
perexchange interaction close to half- filling) , a large effec- 
tive mass enhancement, (to* /m^ = 5 for U = 10|t|), large 
low-T resistivity coefficient, and enhanced Kadowaki- 
Woods ratio. 

Thus our results show that the details of the band 
structure are extremely important for understanding the 
metallic and spin ordered phases of Na^CoC^. In par- 
ticular, a simple one-band Hubbard model on a trian- 
gular lattice is not enough to understand the unconven- 
tional properties of Na^OoC^. Our results however do 
suggest that the Curie- Weiss susceptibility observed in 
the metallic phase of Na a; Co02 results from incoherent 
quasi-localized electrons. Remarkably, the DMFT esti- 
mate of T* = 0.1\t\ s=s 100 — 200 K is in good agreement 
with the Curie- Weiss constant 9 w 150 K experimentally 
measured. 

Many of the thermodynamic and transport properties 
observed in Naa;Co02, for large x, are consistent with the 
DMFT prediction of the transport and magnetic proper- 
ties of the metal very close to the Mott transition. We 
have proposed that Mott physics is relevant to NazCo02 
even when x is large and the system is far away from the 
Mott insulating phase because of experimental observa- 
tion of charge ordering in the Co planes. 

The T dependence of several transport and magnetic 
properties is consistent with the DMFT description of a 
metal close to the Mott transition for T > T* . How- 
ever, at very low temperatures, experimental observa- 
tions on Na x Co02, depart from our DMFT results. For 
instance, the susceptibility deviates significantly from 
the Curie- Weiss behavior being strongly enhanced below 
T ~ 80 K. Also transport and thermodynamic quanti- 
ties are very sensitive to a small applied magnetic field. 
One is tempted to associate the very low temperature 
scale of about 1 K experimentally determined to the 
DMFT coherence scale, T* , which goes to as the sys- 
tem approaches the Mott insulating state. However, such 
a small value of T* implies an extremely large effective 



mass enhancement in clear disagreement with experimen- 
tal values, m* /nib ^3 — 4. We conclude that this very 
small energy scale is due to low energy ferromagnetic 
fluctuations in the Co-planes which are neglected within 
the DMFT. These fluctuations would scatter the strongly 
renormalized quasiparticles which are properly described 
by DMFT. 

We have then used this analysis to help identify 
the simplest relevant model that captures the essen- 
tial physics in NazCo02. We have proposed a model 
which contains the charge ordering phenomena observed 
in the system that we have proposed drive the system 
close to the Mott insulating phase for the large dopings. 
Band structure calculations find that the ai s band can 
be roughly fitted to a tight-binding model with t < 
which roughly describes the large hole pocket around the 
T point. Additional small hole pockets may be present in 
the Fermi surface associated with the e' 2g bands crossing 
the Fermi energy. These bands can be fitted to a t > 
tight-binding band for which ferromagnetic fluctuations 
are likely to exist. Hence, a multiband DMFT treatment 
containing the ai g and e 2g bands as well as charge or- 
dering phenomena should describe both the following as- 
pects: (i) the proximity of the metal to a Mott insulating 
phase leading to renormalized quasiparticles and Curie- 
Weiss behavior for T > T* and (ii) the strong ferromag- 
netic fluctuations appearing at low temperatures. We 
believe that a model including both of these ingredients 
is needed to fully understand the physics of NazCo02. 
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